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Technical  Objectives 

Our  objective  is  to  obtain  a  more  thorough  understanding  of  the  dynamics  of  coupled  air-sea  boundary 
layers  at  relatively  small  spatial  scales,  by  performing  direct  numerical  simulation  (DNS)  and  large-eddy 
simulation  (LE3)  and  targe-wave  simulation  (LWS)  for  both  air  and  ocean  turbulent  flows  with  surface 
waves.  With  extensive  simulation  in  collaboration  with  measurement  we  aim  to  identify  and  assess  the 
key  transport  processes  within  the  atmosphere-ocean  wave  boundary  layer  (WBL).  The  ultimate  goal  of 
this  project  is  to  obtain  the  physical  foundation  for  the  characterization  and  parameterization  of  the 
momentum,  mass  and  heat  transfer  within  the  atmosphere-ocean  WBL. 


Technical  Approach 

For  the  DN5  and  LES/LWS  of  coupled  air  and  ocean  turbulent  flows,  we  have  developed  a  suite  of  high- 
performance,  complementary  computational  methods.  These  include:  (i)  a  boundary  interface  tracking 
method  (B1TM)  for  low  wind  speeds  (<5  m/s);  and  (ii)  an  Eulerian  interface  capturing  method  (EICM)  for 
moderate  to  high  wind  speeds  (>5  m/s),  where  the  waves  steepen/break*  The  numerical  schemes  of  BITM 
are  based  on  boundary- fitted  grids  and  coupled  free-surface  boundary  conditions.  The  EICM  is  based  on 
a  level  set  approach.  These  developments  are  at  the  cutting  edge  of  computational  free-surface 
hydrodynamics.  Transport  of  passive  scalars  in  the  coupled  air-water  flow  system  is  implemented.  Both 
the  BITM  and  EICM  codes  are  optimized  on  parallel  computing  platforms  to  provide  high-resolution 
results  in  a  timely  manner. 

The  BITM  method  solves  the  incompressible  Navier-Stokes  equations  for  both  air  and  water.  Free- 
surface  coupled  boundary  conditions  are  used  at  the  air- water  interface,  with  the  kinematic  boundary 
condition  requiring  that  the  interface  remains  a  material  surface,  and  the  dynamic  boundary  condition 
specifying  a  stress  balance  across  the  interface.  The  transport  of  scalars  is  governed  by  a  convection- 
diffusion  equation.  The  governing  equations  are  discretized  using  a  pseudo-spectral  method  in  the 
horizontal  directions  and  a  finite-difference  scheme  in  the  vertical  direction,  A  second-order  fractional- 
step  scheme  is  used  for  the  time  integration  of  the  flow  field  evolution,  the  transport  of  scalars,  and  the 
motion  of  the  air-water  interface. 

In  EICM,  the  air  and  water  together  are  treated  as  a  system  with  varying  density,  viscosity  and  diffusivity. 
A  continuous  scalar  (the  level  set  function),  which  represents  the  signed  distance  from  the  interface,  is 
used  to  identify  each  fluid.  The  fluid  motions  are  governed  by  the  Navier-Stokes  equations  while  the 
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scalar  is  adverted  with  the  flow  governed  by  a  Lagrangian-invariant  transport  equation.  A  large  wave 
simulation  technique  is  used  to  model  the  effects  of  small  surface  wave  fluctuations  on  large  waves.  The 
governing  equations  are  discretized  on  an  Eulerian  grid  using  a  finite-difference  scheme. 


Work  Completed 

•  Numerical  capabilities 

We  developed  high-performance  DNS  and  LES/LWS  capability  for  coupled  air- water-wave  turbulent 
Hows.  The  numerical  methods  of  RITM  and  EICM  developed  in  this  project  substantially  improve 
the  accuracy  and  physical  meaning  of  the  numerical  results  over  existing  computational  approaches. 
Optimized  codes  for  high-performance  computing  (HPC)  parallel  platforms  have  been  obtained, 

•  Physics  of  air-sea  wave  co  upling 

We  obtained  a  detailed  instantaneous,  three-dimensional  description  of  the  How  field  in  the  coupled 
air- water  turbulent  flow.  The  datasets  obtained  from  our  simulations  provide  us  with  a  physical  basis 
to  gain  useful  insights  into  the  structures  and  dynamics  of  turbulent  flows  in  the  vicinity  of  the  air-sea 
interface.  In  particular,  we  identified  characteristic  vortical  structures  in  the  vicinity  of  the  air-water 
interface  to  elucidate  effects  of  turbulence  coherent  structures  on  transport  processes  in  the  air-sea 
WBL,  Micro  low-level  jets  on  the  air  side  induced  by  water  motion  were  discovered,  which  plays  an 
important  role  in  the  transport  of  passive  scalars  in  the  vicinity  of  the  interface.  We  also  obtained  a 
physical  understanding  of  the  energy  transfer  and  dissipation  process  in  wFater  and  air  during  wave 
breaking  events. 

•  Modeling  of  coupled  air-water  dynamics 

We  quantified  turbulent  kinetic  energy  (TKE)  budget  in  the  vicinity  of  air-sea  interface,  to  establish 
the  physical  basis  of  meso-scale  LES  and  RANS  modeling.  We  also  evaluated  enstrophy  budget  for 
its  parameterization.  Of  significant  importance  is  that  we  quantified  and  completed  model  assessment 
for  surface  renewal  process  that  is  essential  for  air-sea  interfacia  1  transport. 


Results 

In  this  work,  we  have  developed  DNS  and  LES/LWS  capabilities  for  the  coupled  air-ocean-wave  flow 
field,  A  typical  simulation  result  is  plotted  in  Figure  1,  which  shows  the  interactions  among  wind,  swell, 
and  ocean  current.  Such  simulation  provides  important  three-dimensional,  instantaneous  information  on 
the  air-sea  interactions  at  the  sea  surface.  Representative  quantities  including  surface  pressure 
distribution,  transfer  rate  of  momentum  and  energy,  and  vertical  structures  are  shown  in  Figure  L 

Substantial  understanding  of  the  air-sea  coupled  motion  and  the  transport  process  within  the  coupled 
boundary  layers  has  been  obtained  via  the  study  of  a  canonical  problem  of  coupled  air- water  turbulent 
Couette  flow  shown  in  Figure  2,  Through  large-scale,  high-resolution  simulations,  we  are  able  to  obtain 
the  detailed  flow  field  and  scalar  concentration  distributions.  A  typical  example  is  shown  in  the  right  plot 
of  Figure  2,  Based  on  detailed  analyses  of  the  vorticity  field,  characteristics  of  vortices  have  been 
identified.  Vortical  structures  in  the  water  can  be  characterized  as  four  categories:  hairpin  vortices, 
quasi-streamwise  vortices,  interface-attached  single  and  paired  (“LH~shape)  vortices.  The  first  three  are 
illustrated  in  Figure  2.  The  coherent  vortical  structures  identified  in  this  study  are  found  to  play  a 
significant  role  in  the  near-surface  transport  processes. 
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After  the  identification  of  these  instantaneous  structures,  characterizing  and  quantifying  their  effects  on 
near  interface  transport  processes  becomes  a  great  challenge.  In  this  study,  we  employ  a  variable-interval 
space-averaging  (VISA)  technique,  which  is  based  on  the  variable-interval  time-averaging  (VITA) 
method  developed  in  experiment.  For  each  type  of  structures,  0(1 03)  events  are  captured  with  our 
extensive  datasets  of  simulations.  A  typical  example  for  the  quasi-streamwise  vortices  is  shown  in  Figure 
3.  One  important  discovery  from  our  study  is  the  airside  micro  low-level  jets  which  are  induced  by  the 
water  motions  underneath.  The  left  figure  of  Figure  3  shows  that  micro  low-level  jet  is  a  salient  feature  in 
the  region  to  the  upper-left  of  the  vortex.  Also  in  that  region  the  convection  induced  by  the  vortex  makes 
the  thickness  of  the  scalar  boundary  layer  much  thinner,  as  shown  in  the  right  figure  of  Figure  3.  As  a 
result,  the  transfer  of  the  scalar  across  the  interface  is  greatly  enhanced. 

Based  on  the  extensive  simulation  datasets,  we  obtained  important  statistics  for  the  modeling  and 
parameterization  of  atmosphere-ocean  coupled  boundary  layer.  In  this  study  we  are  able  to  complete  the 
TKE  budget  in  atmosphere-ocean  WBL,  which  is  governed  by  a  number  of  processes.  The  TICE  is 
produced  by  the  interaction  between  Reynolds  stress  and  the  shear  in  the  mean  flow.  Meanwhile,  it  is 
dissipated  by  viscosity.  In  addition  to  these  source  and  sink  processes,  TKE  is  transported  among 
different  flow  regions  by  pressure  fluctuations  and  velocity  fluctuations.  At  regions  where  the  TKE 
varies  abruptly,  viscous  diffusion  also  contributes  to  the  TKE,  Figure  4  shows  that  the  TKE  production  is 
large  near  the  interface.  The  sharp  decrease  at  the  interface  is  caused  by  the  reduction  of  Reynolds  stress 
there.  Viscous  diffusion  is  only  significant  very  close  to  the  interface.  Transport  due  to  turbulence 
velocity  fluctuations  transports  TKE  from  the  bulk  region  of  the  air  to  the  near-surface  region.  On  the 
waterside,  turbulence  transport  removes  part  of  TKE  from  the  near-surface  region  and  put  it  at  the  deep 
region.  Dissipation  increases  towards  the  interface  in  the  air  and  reaches  a  maximum  at  the  interface, 
behaving  like  near  a  solid  wall.  On  the  waterside,  as  the  interface  is  approached,  dissipation  decreases 
first  and  then  increases.  The  results  obtained  here  will  allow  us  to  develop  advanced  closure  models 
needed,  for  example,  in  larger  scale  and  coarser  grid  simulations  such  as  those  used  in  meso-scale  LES 
and  RANS  computations. 

We  also  obtained  valuable  statistics  of  surface  renewal,  which  is  the  most  important  interfacial  transport 
process.  Of  significant  importance  is  that  we  develop  a  novel  Lagrangian  and  Eulerian  approach  for 
quantifying  ages  of  fluid  surface  elements.  Figure  5  shows  the  probability  density  obtained  from  our 
direct  numerical  simulation.  Various  theoretical  results,  including  exponential.  Gamma,  normal,  and 
lognormal  distributions  are  examined.  It  is  found  that  the  distribution  of  surface  age  is  not  Gaassian,  The 
exponential  distribution  based  on  the  assumption  of  random  surface  replacement,  which  is  widely  used  in 
the  literature,  does  not  capture  the  entire  process  of  the  surface  renewal.  For  young  surface  elements,  the 
lognormal  is  a  good  fit,  while  for  old  elements  the  Gamma  curve  works  better.  Considering  that  the 
majority  of  surface  flux  occurs  right  after  surface  renewal,  we  conclude  that  the  surface  age  distribution 
can  be  best  described  by  the  lognormal.  This  result  establishes  a  physical  basis  for  the  parameterization 
of  scalar  transport  at  the  air-sea  interface. 

Impact/Application 


This  study  aims  to  obtain  a  fundamental  understanding  of  the  air-sea  wave  coupling  dynamics  at  small 
scales  at  low  wind  speeds.  Our  work  is  a  small  yet  essential  part  of  an  overall  coordinated  effort 
involving  field  experimentalists,  air-sea  modelers,  and  physical  oceanographers  to  obtain  improved 
physics-based  parameterizations  for  air-sea  interactions.  Our  numerical  simulations  provide  detailed 
descriptions  of  the  air-sea-wave  boundary  layer  at  small  scales,  and  a  physical  basis  for  the  modeling  and 
parameterization  of  transport  process  within  the  atmosphere-ocean  wave  boundary  layer.  The  simulations 
also  provide  comparison  and  cross-validation  with  field  measurements.  Finally,  our  numerical  framework 
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can  be  used  as  a  powerful  tool  to  help  the  interpretation  and  syntheses  of  field  data,  and  parameterization 
of  WBL  transport  process  for  the  coupled  air-ocean  WBL  modeling* 
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Figure  1.  Interaction  among  wind,  swell,  and  underwater  current.  Plotted  on  the  left  are  surface  wave 
profile  and  contours  of  air  and  water  velocity  in  the  direction  of  wave  propagation*  The  contours  on 
the  right  are  surface  distributions  of  vertical  vorticity,  pressure,  momentum  and  energy  inputs  from  the 

atr  to  water,  respectively. 
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Figure  2.  Left:  schematics  of  air-water  coupled  Couette  flow .  Right:  instantaneous  vortical  structures 
in  the  air  and  water  turbulent  flow  field ,  and  instantaneous  scalar  concentration  contours. 
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Figure  3.  Quasi-stream  wise  vortex  on  the  waterside  and  the  induced  mirror  vortex  and  micro  low-level 
jet  on  the  air  side,  and  their  effects  on  passive  scalar  transport  near  the  interface .  Shown  in  the  left 
figure  are  contours  of  vorticity  component  and  fluctuation  velocity  vectors .  Shown  in  the  right  figure 

are  contours  of  scalar  concentration  fluctuations. 


Figure  4.  Profiles  of  turbulence  kinetic  energy  budget  terms  in  the  coupled  air-water  boundary  layer. 


6 


Figure  5.  Probability1  distribution  of  surface  element  age  with  respect  to  the  surface  renewal  process. 
Histogram  is  obtained  from  simulation  data .  Curve  fittings  are  exponential ,  Gamma,  lognormal ,  and 

normal t  respectively. 
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